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Phonon dispersion in a two-dimensional electron-lattice system described by a two- 
dimensional square-lattice version of Su-SchriefFer-Heeger's model and having the half- 
filled electronic band is studied theoretically at temperatures higher than the mean field 
critical temperature of the Peierls transition. When the temperature is lowered from 
the higher region down to the critical one, softening of multi phonon modes which have 
wave vectors equal to the nesting vector Q = {n/a, ir/a) with a the lattice constant or 
parallel to Q is observed. Although both of the transverse and longitudinal modes are 
softened at the critical temperature in the case of the wave vector equal to Q, only the 
transverse modes are softened for other wave vectors parallel to Q. This behavior is 
consistent with the Peierls distortions at lower temperatures. 

KEYWORDS: two-dimensional electron-lattice system, phonon softening, Peierls transition, longitunidal 
mode, transverse mode, linear mode analysis, SSH model 



1. Introduction 

The Peierls transition'^^ which is caused by the competition between the energy loss due to 
lattice distortions and the gain due to the formation of a gap at the Fermi level in the electronic 
energy spectrum is one of peculiar properties of one-dimensional electron-lattice systems. It is 
also well known that this transition can occur even in higher dimensions if the structure of the 
Fermi surface allows strong nesting by a single nesting vector. A two-dimensional square lattice 
system with a half-filled electronic band described by a tight-binding model is a typical example 
where the Fermi surface (more properly Fermi line) can be completely nested by the wave vector 
Q = (7r/a,7r/a) with a the lattice constant. Early studies on the issue of two-dimensional Peierls 
transition in the half-filled systems suggested that the lowest energy phase would be the state with 
lattice distortions whose Fourier components involved only Q-modes."^^ The relation between the 
electron-electron interactions and the two-dimensional Peierls instability has been discussed rather 
intensively in connection to the mechanism of high Tc superconductivity. ■'' Nevertheless 

the investigation of the true ground state of the pure electron-lattice system without electron 
interactions was not clearly understood until recently. In recent works by Hamano and one of the 
present authors,'^' "^^ it has been shown that the lowest energy state of a two-dimensional electron- 
lattice system described by a two-dimensional version of Su, SchriefFer and Heeger's (SSH) model 
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Hamiltonian') involves not only the distortions with Q-vector but also those with various wave 
vectors parallel to Q, the latter contributing to the formation of the electronic energy gap through 
a second order process (see Fig. ??), and that there are an infinite number of non-equivalent lattice 
distortion patterns which indicates an infinite degeneracy of the lowest energy states. It has also 
been argued'^ that those behaviors are maintained at finite temperatures and that the amplitudes 
of lattice distortions with different wave vectors, all being parallel to Q, vanish all together at a 
critical temperature. 




— TT 



Fig. 1. The situation of the Fermi-hne nesting in the half-filled two-dimensional square-lattice system. The thick-line 

square represents the Fermi line, and a solid line is the nesting vector connecting two different parts of the Fermi 
line. The dotted lines describe the nesting due to the second order perturbation, where two different parts of the 
Fermi line are connected through a second order process. 

It is also known^^ that the Peierls transition in one-dimensional systems is accompanied by the 
softening of the 2k^ phonon mode at the transition temperature, where k-p represents the Fermi 
wave number and 2k-p corresponds to Q in the two-dimensional square lattice systems with a half- 
filled electronic band. As has been clarified by previous papers, phonon modes relevant to the 
Peierls transition in the two-dimensional square lattice system described by the SSH model involve 
not only the Q modes but also many other modes with wave vectors parallel to Q. Furthermore, 
as discussed in refs. ? and ?, the lowest energy states are not unique. In this sense it will be 
worthwhile to investigate the phonon dispersion at finite temperatures in this two-dimensional 
system. Even at temperatures higher than the transition temperature, we may expect a phonon 
softening phenomenon which is much different from that in one-dimension. In the present paper, 
we consider the phonon dispersion of the two-dimensional square lattice system with a half-filled 
electronic band only in the temperature region higher than the critical one, and study which types 
of phonon modes are softened when the temperature approaches the critical one from above. The 
dispersion at lower temperatures cannot be treated in a simple way because of the degeneracy of 
the lowest energy states, and therefore it will be discussed in a forthcoming separate paper. 

Since the formation of the energy gap at the Fermi level leads a metallic system into an 
insulator, it is quite important to understand the mechanism of the Peierls transition. The present 
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work will shed light on the understanding of the Peierls instability in two-dimensional electron- 
phonon systems. 

The paper is organized as follows; in the following section the model for and the method of the 
calculation are described. In Sec. 3, the results of the calculations are shown. The last section is 
devoted to summary and discussion. 



Model and Formulation 



In the following calculations we use a two-dimensional version of the SSH model^^ which 
was originally introduced in the studies of one-dimensional electron-lattice systems and where the 
electron-lattice coupling is introduced through the dependence of the electronic nearest-neighbor 
transfer integral on the hopping distance. Since we are treating a two-dimensional square lattice, 
the explicit form of the model Hamiltonian is written as follows, 



H = 



[to - a{ux{i + l,i) - Ux{i,j))] 
+ [to - a{uy{i,j + 1) - Uy{i,j))] 

+ {uy{i,j + 1) - Uy{i,j)) 
M 



+- 



(1) 



where the field operators Ci^j^s and cj^-^ annihilate and create an electron with spin s at the site 
respectively, and to is the transfer integral for the equidistant lattice, a the electron-lattice 
coupling constant, u{i,j) the lattice displacement vector whose x- and y-components are denoted 
as Ux{i,j) and Uy{i,j), respectively, K the force constant describing ionic coupling strength in the 
lattice system, M the mass of an ion unit at a site. Through out this paper, we assume the periodic 
boundary conditions (PBC) for both directions. 

The equations of motion for the lattice system can be derived on the assumption that the 
ionic mass M is much larger than the effective mass of an electron and by introducing the Born- 
Oppenheimer approximation;') in the case of the x-components of the displacement vectors, they 
are explicitly written as follows. 
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+K{ux{i + 1, j) - 2ux{i,j) + ux{i - 1, j)), (2) 

and similar expressions for the y components can be written down, which will be omitted here 
because they can be easily speculated from the above ones. Assuming the electronic system is in 
the equilibrium state for the given lattice configuration at any moment, the above equations of 
motion can be rewritten as follows, 

Mux{i,j) = 2a^/(£^)U4i + l,j>^(z,j) 

+K(ux{i + - 2ux{i,j) + Ux{i - l,i)), (3) 

where f{£u) represents the Fermi distribution function for the cigenenergy corresponding to 
the single particle eigenstate (f)^, which is obtained by solving the eigenvalue problem for a given 
configuration of the lattice displacements { u{i,j)}. Because of the electron-hole symmetry of the 
present system, the electronic chemical potential is fixed to be zero at any temperature. Again the 
equations for the y-component are omitted. 

The stable static lattice displacements { u^{i, j) } are obtained by solving a set of self-consistent 
equations which are nothing but those given by setting the right-hand side of eq. (??) to be zero. The 
linear modes around this static solution, i.e. the phonon modes, can be calculated by linearizing 
the equations of motion for Ux{i,j) and Uy(i,j) with respect to the deviations from the static 
configurations. In deriving equations for the linear modes, we should take into account the effects 
of those deviations on the electronic eigenfunctionsr^ 

u{i,j,t) = vP{i,j) + Su{i,j,t) (4) 
^v{i,3;t) = 4>l{i,j) + S(l)^{i,j;t), (5) 

where we have introduced the time variable t explicitly and the electronic eigenfunctions corre- 
sponding to the static displacements are denoted as (p^{i,j). The time dependence of the electronic 
wave functions is implicit one through the slow time dependence of the lattice displacements. The 
deviations of the wave functions 5(j)i,{i,j\ t) are calculated within the first order perturbation theory 
in the form, 

dMhj;t)= Yl (6) 
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with 



= 



5ux{n + l,m,t) — Sux{n, m, t) 



x^<(n+l,™«(„,™)^ 

+ ^5uy{n, m + l,t) — 5uy{n, m, t)^ 

x(<Pnin,m + l)(f)l{n,m) 



(7) 



In order to determine the normal modes, we postulate that the time dependence of the devia- 
tions of the lattice displacements is given by 



du{i,j,t) = {5ux{i,j,t),duy{i,j,t)) = 6u{i,j;u)e"^* . 
Then the equation for du{i,j;Lu) can be expressed in the following form, 

Mu'^6u{i, j;uj) = W(^, j', m, n)Su{m, n; uj), 

m,n 

where the 2x2 matrices W{i,j;m,n) are defined as follows. 



(8) 



(9) 



bl{r)\^^l{r + ea) -<!>>- ea) 
+4>l{r)(4>l{r + ea)-ct>l{r-ea 



l{R)(4>l{R + e,)-^l{R-e^) 



+4>l{R){cl,l{R + e,) 



(10) 



Here a and h stand for x or y, and for simplicity we expressed the lattice sites as r[= and 
R[= (m, n)\ and introduced unit vectors Cx = (1, 0) and By = (0, 1). The eigenfrequency lo and the 
corresponding cigcnfunction Su{r; uj) for the normal modes are determined by solving the eigenvalue 
problem eq. (??)• 

The above-mentioned formulation can be used at any temperature irrespective of the presence 
or absence of finite static lattice distortions. In the present paper we consider only the high temper- 
ature region where there is no static lattice distortion. In this situation, the normal modes equations 
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can be reduced into a much simpler form. If there is no distortion, the electronic eigenenergies and 
eigenf unctions are given by those of the plane waves; 



£fe = -2to[cos(A;j;) +cos(A;j/)], (11) 
= (12) 

where k = {kx,ky) is the wave vector characterizing the electronic cigcnstatc, and the system size 
is assumed to be L x L. We use a length unit where the lattice constant is equal to unity. It can 
be easily confirmed that the phonon normal modes arc also expressed in the plane wave form in 
the present situation, i.e. 6u{r,Lo) = g{q,Lo)e^'^'''^ . Thus the linear modes equations are rewritten 
in the following form, 

uj'^g{q,uj) =U{q)g{q,Lu), (13) 



where the elements of the 2x2 matrix U{q) are given by 

k 



Ua,b{<l) = -TJJ2 y] 

ML^ - Sk- £k+q 



X |^sin(A:a + qa) - sin ka j \ sm{kh + %) - sin k^ 
K 

+ cos qa)5a,h- (14) 

Here again a and h stand for x oi y. It is easily confirmed that lA{q) is a symmetric matrix. The 
eigenfrequency lo does not appear in the matrix U{q) because of the Born-Oppenheimer approx- 
imation. It is straightforward to calculate the eigenvalues and eigenvectors from eq. (??). The 
results are discussed in the following section. It should be noted that the calculation of the phonon 
dispersion in the low temperature region, where distortions involving various wave numbers parallel 
to the nesting vector Q are appearing, is not so simple as in the present situation without any static 
distortion. 



3. Results of Calculations 

As stated in the previous section, the frequencies of the phonon modes in the temperature 
region higher than the Peierls transition temperature are obtained from the eigenvalue problem 
for a 2x2 matrix W, the elements of which depend only on the wave vector q and the temperature 
included in the electronic distribution function (Fermi distribution). The result is written as 

= ^ [U:,,X + Uy,y ± ^{U,,,-Uy,yf+mly] (15) 

When Qx = Qy = 0, all the elements of the matrix U vanishes irrespective of the temperature 
and consequently we have two degenerate uniform modes with zero frequency; they correspond 
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to the long- wave-length limit of two acoustic modes. For general values of Qx and qy, the matrix 
elements are calculated only numerically. 

As has been shown in the previous works,^' ""^ the phonon modes relevant to the Peierls transition 
are those with wave vectors parallel to Q including Q itself. Therefore, let us first discuss the 
temperature dependence of several phonon frequencies for the case with Qx = Qy = Q- In this 
case, the phonon modes are longitudinal {g{q) || g, i.e. gx = Qy) or transverse {g{q) -L q, i.e. 
gx = —gy)V In Fig. ?? the eigenvalue a;^ is plotted as a function of the temperature T for the 
transverse (a) and longitudinal (b) modes with q = tt (solid curve), 7r/2 (dashed curve) and 7r/4 
(dotted curve) . The product of the temperature T and the Boltzmann constant is scaled by to and 
a;^ by K/M. The dimensionless coupling constant A = a^/toK is assumed to be 0.32 throughout 
the paper. In evaluating the fe-sums included in the definition of the matrix elements Ua,b, we have 
set the system size L to be 64; this choice of the system size is to get a clear presentation of the 
wave vector dependence of the dispersion (see Fig. ??). It has been confirmed that the results 
shown in Fig. ?? do not change even if larger values of L are used. 

(b) 
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Fig. 2. The temperature dependence of the eigenvalue uP' for transverse (a) and longitudinal (b) modes with fixed 
values of the wave number q {= qx = Qy)', the values of q axe given in the graph. In the inset of (a) all eigenvalues 
for the transverse modes with Qx = qy = q (including q = tt ) scaled by corresponding values for free phonons u>o (q) 
are plotted; they fall into a single curve. The temperature [ksT with fee the Boltzmann constant) is scaled by to 
and u!^ by K/M. The dimensionless coupling constant is A = 0.32 and the system size L is chosen to be 64. 



In Fig. ??(a), there are three examples of temperature dependence of the eigenfrequencies for 
the transverse modes with wavevectors qx = Qy = Q are plotted, it should be noted that all the 
eigenvalues cross at the same temperature (fceTc ~ 0.35). We have numerically confirmed that the 
eigenvalues u>'^{q,T) for the transverse modes with Qx = Qy = Q scaled by uiq [= {K/M){1 — cosqr)] 
the square frequency of the corresponding free phonon depend on the temperature T but not on q 
as shown in the inset of Fig.2 (a). It will be clear that all the frequencies of the transverse modes 
with wavevectors Q and the parallel to it vanish at the same temperature. This is consistent with 
the structure of the lower temperature Peierls phase.'' The critical temperature Tc is the same 
with that estimated from the vanishing of the order parameters in the low temperature regime.''^ 



8 Shutaro Chiba and Yoshiyuki Ono 

In the case of the longitudinal modes, only the mode with the nesting vector (i.e. q = tt) vanishes 
at the critical temperature as seen in Fig. ??(b). When the temperature is reduced below the 
critical temperature, the values of for all the modes in Fig. ??(a) and that for the nesting vector 
mode in Fig. ??(b) become negative, indicating these modes are unstable in the lower temperature 
regime and suggesting a phase transition. 

In order to see the softening of the phonon modes with wave vectors parallel to the nesting 
vector at the critical temperature, we have calculated the wave vector dependence of the phonon 
dispersion at Tp. The results are summarized in Fig. ??, where in (a) the dispersion of the lower 
branch [see eq. (??)] which corresponds to the transverse modes for some special directions of the 
wave vector is shown and in (b) the upper branch connected to the longitudinal modes for some 
special directions. 

For the sake of simplicity we call the lower branch [Fig. ??(a)] the transverse modes and the 
upper branch [Fig. ??(b)] the longitudinal modes. In the case of the transverse modes, all those 
with wave vectors parallel to the nesting vector Q, including the Q-mode itself, are found to have 
zero frequency at the critical temperature, while in the case of the longitudinal modes, only the 
frequency of the Q-mode vanishes at the critical temperature. 

4. Summary and Discussion 

In the present work, the temperature dependence of the phonon dispersion in the high temper- 
ature regime (i.e. at temperatures higher than the Peierls transition temperature) has been studied 
for a two-dimensional electron-lattice system with a half-filled electronic band, which is described 
by the SSH-type model extended to the case of a two-dimensional square lattice. In this system, it 
was shown in the previous papers ^' that the lowest energy state in the low temperature regime 
is the Peierls state with multi-modes distortions where the Fourier components of the distortions 
involve not only the nesting vector modes but also many other modes having wave vectors parallel 
to the nesting vector Q = (vr, vr). In accordance with this low temperature behavior of the system, 
we have found the softening of phonon modes with wave vectors Q and with those parallel to it 
as the temperature is lowered from the high temperature region towards the critical temperature. 
As for the polarization of the softened modes, both of the longitudinal and transverse modes are 
softened in the case of the Q-vector, while only the transverse mode has zero frequency in the case 
of wave vectors parallel to but not equal to Q at the critical temperature. 

Although the zero temperature ground state of the system has an infinite number of degen- 
eracy,^-' depending on different amplitude distributions of the Fourier components of the static 
Peierls distortions, we have confirmed that the polarization of distortion components correspond- 
ing to non-Q-modcs (i.e. those with wave vectors parallel but not equal to Q) is transverse and 
that the polarization of the Q-mode is longitudinal, in all cases. This fact is maintained even at 
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finite temperatures (below the critical temperature). In spite of the softening of the transverse 
Q-mode, there is no transverse component of the Q-mode in the ground state Peierls distortions. 
On the other hand, the asymmetrically dimerized Peierls state, which was suggested by Tang and 
Hirsch') but has a higher energy than the multi-modes Peierls state,') involves the longitudinal 
and transverse components of the Q-mode with an equal amplitude in the Fourier expansion of the 
distortions. 

Similarly, although all the transverse modes with wave vectors parallel to Q are softened, not 
all the components appear in the Fourier expansion of the multi-modes Peierls distortions. The 
simplest pattern in the degenerate ground states will be the one consisting of the longitudinal 
Q-mode and the transverse Q/2-mode.'') These facts indicate that, when many phonon modes 
are softened, it is not necessary for all the modes to show condensation in order to realize the 
lowest energy Peierls state. At the moment, there is no guiding principle to determine which of the 
infinitely degenerate low temperature states is chosen when the temperature is lowered from the 
higher region across the critical temperature. 

It will also be worthwhile to study the phonon dispersion in the Peierls distorted phase, which 
is to be treated in the future work. Furthermore the effect of electron-electron interactions should 
be considered in order to discuss the Peierls transition in real systems, which is left for future 
studies. 
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Fig. 3. Phonon dispersions at Tc. (a) the dispersion related to "transverse" modes and (b) that related to "longi- 
tudinal" modes. The eigenvalue uP' is scaled by KIM. The dimensionless coupling constant A is 0.32, the system 
size L being 64. 



